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Abstract 

It is known that there exist 32 triplets of circles such that each circle 
is tangent to the other two circles and to two of the sides of the triangle 
or their extensions. We provide formulae to obtain the radii of the circles 
for each of the 32 triplets from the side lengths of the reference triangle 
by means of trigonometric or hyperbolic functions. 

1 Introduction 

The configuration of three circles inside a triangle such that each circle is tangent 
to the other two circles and to two of the sides of the triangle has been studied 
for more than two centuries. Today, such three circles are called the Malfatti 
circles of the triangle. 




Sometime before 1773, Naonobu Ajima (17327-1798), who was a samurai, 
or a member of the military class in old Japan, found a method to calculate 
the diameters of the Malfatti circles from the side lengths of an arbitrary trian- 
gle. The method was called Nanzan-shi san-sha naiyo san-en jutsu ("Nanzan's 
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method on a triangle that includes three circles" , as Nanzan is a pen name of 
Ajima's) or San-sha san-en jutsu in short. A brief description of the method 
is found in 1^ I ^14]. A detailed description of the method including a proof 
is found in [T]. Unfortunately, Ajima's method as well as any other results by 
Japanese mathematicians in those days was inaccessible from outside Japan un- 
til the Edo shogunate, the former government of Japan (1603-1868), abandoned 
the isolation policy in 1854. 

In 1803, an Italian mathematician Gianfrancesco Malfatti (1731-1807) [TU] 
gave a construction to draw the Malfatti circles for an arbitrary triangle. Despite 
Malfatti's unawareness of Ajima's works, Malfatti's construction is considered 
identical in many parts to Ajima's method. 

In 1852, Schellbach [T^[T3] gave a set of formulae to obtain the distances 
between the vertices and the tangent points of the circles on the sides from 
the side lengths of an arbitrary triangle by using trigonometric functions. The 
same formulae with a proof essentially identical to Schellbach's are described in 
English in §30] [7]. 

In 1895, Derousseau [3] generalized the Malfatti circles by removing the 
condition that the three circles are inside the triangle. Derousseau proved that 
there exist 32 triplets of circles such that each circle is tangent to the other 
two circles and to two of the sides of the reference triangle or their extensions. 
Some alternative proofs of the existence are known. In 1904, Pampuch [11] gave 
another proof. In 1930, Lob and Richmond [5] gave yet another proof. 

In this article, we provide formulae to obtain the radii of the circles for each 
of the 32 triplets from the side lengths by means of trigonometric or hyperbolic 
functions. In other words, we provide Schellbach-style formulae for all of the 
Derousseau-Pampuch generalizations. 

2 Notation 

Throughout this article, we use the following notation. 

For a triangle ABC, let a, b, c denote the lengths of the sides BC, CA, AB, 
s the semiperimeter, r the inradius, and ta, tb, re the exradii as usual. 

Let the incircle is tangent to the side BC at D, to the side CA at E, and to 
the side AB at F . Let the excircle corresponding to the vertex A is tangent to 
the side BC at Da, to the extension of the side AC at Ea, and to the extension 
of the side AB at Fa ■ 

Suppose the circle A'{ri) is tangent to the line CA at Ei and to the line 
AB at Fi, the circle B'{r2) is tangent to the line AB at F2 and to the line BC 
at £'2, the circle C'{r^) is tangent to the line BC at D3 and to the line CA at 
£3, and the three circles are tangent to one another. Suppose the nine tangent 
points are distinct. 

3 Classification 

Since the center A' does not locate neither on the line AB nor on the line AC , 
it locates inside ZCAB, inside ZCAB, inside ZCAB or inside ZCAB where 
an overline indicates that the angle has, as one of its sides, the opposite ray 
instead of the ray including the triangle side. For example, ZCAB denotes the 




Figure 2: 



angle with the ray AC and the ray opposite to the ray B . And ZCAB denotes 
the vertical angle of ZCAB. Analogously, the center B' locates inside ZABC, 
inside ZABC, inside ZABC or inside ZABC and the center C locates inside 
ZBCA, inside ZBCA, inside ZBCA or inside ZBCA. 

If the circles A'{ri), B'{r2), C"(r3) lie in Z\i, A2, A3, respectively, then 
/ii n Z\2 ¥" 0: ^1 n Z\3 ^ 0, and Z\2 n Z\3 ^ since the three circles are 
tangent to one another. Thus, for locations of the three centers A', B' , C , only 
7 out of the 64 cases are consistent to the condition that the circles are tangent 
to one another. The following are the consistent cases. 





A' is inside 


B' is inside 


C is inside 


Case 1 


ZCAB 


ZABC 


ZBCA 


Case 2 


ZCAB 


ZABC 


ZBCl 


Case 3 


Z'CAB 


ZAB'C 


Z'BCA 


Case 4 


ZCAB 


ZABC 


ZBCA 


Case 5 


ZCAB 


ZABC 


ZBCA 


Case 6 


Z'CAB 


ZAB'C 


ZBCA 


Case 7 


ZCAB 


ZlBC 


Z'BCl 



4 Solution 

4.1 Case 1 

In Case 1, the following three conditions hold. 



BD2 + D3C + D2D3 ^ 


= BD + DC or BD2 + D3C - 


- D2D3 


= BD + DC 


AEi + EsC + EiEs-^ 


^ AE + EC or AEi + E3C - 


- E1E3 = 


= AE + EC, 


AFi + F2B + F1F2 


^ AF + FB or AFi + F2B - 


-FiF2 = 


-- AF + FB. 



By expressing the lengths by the radii and the angle sizes, we obtain from the 
first disjunction that 



B 



C 



B 



C 



r2 cot y + ^3 cot — + 2^/r^ = r cot — + r cot — 



or 



B 



C 



B 



C 



r2 cot — + rs cot — - 2^/r^ = r cot — + r cot — , 



we obtain from the second disjunction that 



A 



C 



A 



C 



ri cot — - + r3 cot — + 2^rir^ = r cot tt + r cot ■ 



or 



A C „ , A 

ri cot h rs cot zVr-ira — r cot \- r cot 

2 2^ 2 



C_ 
■2' 



and we obtain fr-om the second disjunction that 



A 



B 



A 



B 



n cot -^ + ^2 cot — + 2^rir2 = r cot — + r cot — 



or 



B 



B 



ri cot h r2 cot 2^r\r2 = r cot h r cot — . 



(1) 

(2) 

(3) 
(4) 

(5) 
(6) 



Define l^ m, nhy 

A 
I = cot-, 

Define u, v, w, x, y, z by 



B 

m — cot — , 
2 



Then we have 



r 


if dH) holds, 


< 


if @ holds. 


r 


if holds. 


< 


if dH) holds. 


( y/rir2 


if © holds, 


1 V^l'^2 


if ra VinlHs 



J"! 



r 
V - 



C 
n — cot — . 
2 



r 



r3 



' my + nz + 2u = m + n, 
Ix + nz + 2v ^ I + n, 
Ix + my + 2w = I + m,, 
xy — w^, 
xz = v^ , 
yz — u . 



(7) 



(8) 



For any triangle ABC, if I, m, n are defined by (O, then Irrm = I + m + n 
holds. On the other hand, if positive reals I, m, n satisfy Imn — I + m + n, 
then there exists a triangle ABC that satisfies ([T]). Thus, Case 1 can be reduced 
into solving the system of equations ^ for u, v, w, x, y, z with positive real 
parameters I, m, n under the restriction Imn — I + m + n. 

As we will show in Appendix |3 the system of equations has the following 8 
solutions. 



w 



VP 


+ 1- 


/ + 1 




2 


' 


y/m 


2 + 1 


-m + 1 




2 




y/n^ 


+ 1- 


-n+1 



I + m + n - 


- 1 + \/P + 1 - Vm^ + 1 - Vn'^ + 1 




21 


I + m + n - 


-l-Vl'^ + l + Vm^ + 1 - Vn^ + 1 




2m 


I + m + n - 


-1-Vl^ + l- Vm^ + 1 + Vn^ + 1 



2n 



VP +1-1+1 






2 




\Jm?- + 1 + 771 


- 1 




2 




\rn 


2 + 1 + 77- 


1 



/ + 777 + 77 - 1 + VWTT + \Jm^ + 1 + \/«2 + 1 



v = 





- 1- 


21 


I + m. + n - 


-VP + l- Vm^ + 1 + Vn2 + 1 




- 1- 


2777 


/ + 777 + 77 - 


-Vl^ + l + Vm2 + 1 - Vn^ + 1 



2n 





V^2 + l+/-l 




2 




y/m'^ + 1 — 777+1 




2 



Vr72+T+77- 1 



/ + 777 + 77 - 1 - VWTT - \/m2 + 1 + V?l2 + 1 



y = 





21 


I + m. + n - 


- 1 + V^2 + 1 + ^jn^ + 1 + ^„2 _^ 1 




2777 


I + m, + n - 


-l + ^/P + l- Vm2 + 1 - V«2 + 1 



277 



(9) 



(10) 



(11) 



VWTT + i-1 



\Jm?- + 1 + m — 1 



w 



\/'n? + 1 — n + 1 



/ + m + n- 1 - VF + T + Vto2+T - Vn2+T 





21 


I + m + n - 


-l + Vl^ + l- Vm^ + 1 - Vn2 + 1 




2™ 


I + m + n - 


-l + Vl^ + l + Vm^ + 1 + Vn^ + 1 



2n 



v;2 


+l+;+l 


2 


Vn^ 


2 + 1 + m+l 


2 


V"-^ 


+ l+n+ 1 



y = 



/ + 


m + 


n+l-VF 


+ 1 + Vm"^ 


+ 1 + Vn2 


+ 1 




m + 


n + l + VF 


21 






; + 


+ 1 - \/?Tl2 


+ 1 + Vn2 


+ 1 




m + 


n + l + VF 


2m 






/ + 


+ 1 + Vto2 


+ 1 - Vn2 


+ 1 



2n 



VI 


^+l+?+l 




2 


\/m?- 


+ 1 -m- 1 


2 


Vn2 


+l-n-l 



y = 



/ + 


m + 


n + l-VF 


+ 1- 


Vm2 


+ 1- 


-Vn2 


+ 1 




m + 




21 










/ + 


n + l + VF 


+ 1 + Vto^ 


+ 1- 


-Vn2 


+ 1 




?Ti + 


n + l + Vl2 


2m 










/ + 


+ 1- 


Vm2 


+ 1 + Vn2 


+ 1 



2n 



(12) 



(13) 



(14) 



VpTT-1-1 

2 ' 

Vm^ + 1 + m + 1 



w 



\/v? + 1 — n — 1 



/ + 


m + 


n + l + Vl2 


+ 1 + Vm.^ 


+ 1- 


-V^ 


+ 1 




?n + 


n+1- 


-VF 


2? 








/ + 


+ 1 - v™^ 


+ 1- 


-Vn^ 


+ 1 




m + 


n + 1 - 


-x/I^ 


2™ 








/ + 


+ 1 + \/m^ 


+ 1 + Vn^ 


+ 1 



2n 



V^2+T-;-i 



Vm^ 


+ 1 -m 


-1 




2 


; 


Vn 


2 + 1 + n 


+ 1 



Z + m + n + 1 + \/7M~T- Vrn.2 + 1 + V"^ + 1 



y = 





2? 


/ + m + n + 1 - 


- V^2 + 1 + ^7^2 + 1 + Vn^ + 1 




2m 


/ + m + n + 1 - 


- \//2 + 1 - Vm2 + 1 - Vn2 + 1 



2n 



(15) 



(16) 



Define a,/3 


7G(0 


7r/2) and a by 












. 2 

sm a = 


a 

1 
s 


• 2 n 

sm /3 = 


b 

s' 


sin^-) 


c 


a = 


a + (5 + 
2 


7 


The fourth 


equation in ^ 
r(^l + m + n - 


corresponds to 


a value of 


ri as 


follows. 






l + x//2 


+ 1- 


Vm2 + 1 - 


-Vn^ 


+ 1) 





Since 



2/ 



(s — a){s — 6)(s — c) 



-6 



r rA 



r tb 



re 



it holds that 



r{l + m + n-l + VPTT - \/m? + 1 - \/n^ + l) 



21 




(s — a)(s — \))(s — c) Us — a)bc a{s — b)c ab{s — c) 



— (1 — cos a cos /3 cos 7 + cos a sin /3 sin 7 — sin a cos /3 sin 7 — sin a sin /3 cos 7) 
r/i(l - cos(;3 + 7 - a)) 



= r^ sin (cr — a). 

By making similar calculations on every last three equations in ^, ([TU)) . 
(fTTj) . (fT2]) . ([T3)) . ([T4|. (fT5)) and dHl), we obtain the following respective solutions 
in Case 1. 



ri = rA sin^(cr — a), 

'"2 = rs sin^ (cr - /3), 

rs = rcsin^(cr - 7). 

ri = r/i sin cr, 

r2 = rs sin^(cr- 7), 

rs = rcsin^(cr - ^). 

ri = ?:4sin^(cr-7), 

'^2 = rs sin^ cr, 

''a = i"c sin^(cr — a). 

ri = rAsin^(cr- /3), 

r2 = rBsin^(cr- a), 

r3 = rp sin cr. 
' ri = rA cos^ (cr — a), 

r-2 = rBCOs^(cr- /3), 
Kr3 = rccos^ia - 7). 

ri — rA cos^ cr, 

?'2 = rs COS^(cr- 7), 

rg = re cos^(cr- /3). 

ri = rA cos^(ct — 7), 

?'2 = tb cos^ ct, 

r3 — re cos^(cr — a). 
' ri = rACOS^(cr-^), 

7'2 = rB cos^(ct — a), 
,r3 = rc cos^ cr. 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



4.2 Cases 2 & 3 

In Case 2, we have that the fohowing three disjunctions of equations hold. 

B£»2 + D3C + D2D3 = BDa + DaC or BD2 + D3C - D2D3 = BDa + Da C, 
AEi - CE3 + E1E3 = AEa - CEa or AEi - CE3 - E1E3 = AEa - CEa, 
AFi - BF2 + FiF2= AFa - BFa or AFi - BF2 - F1F2 = AFa - BFa. 

By expressing the lengths by the radii and the angle sizes, we obtain from the 
first disjunction that 



B C „ , B C 

r2 tan — + rs tan — + 2V?'2?'3 = ^4 tan "^ + ^4 tan — 



or 



B 



C 



^ „ , , B C 

r2 tan — + rs tan — - 2yjr2rz = rA tan — + rA tan — , 



we obtain from the second disjunction that 



A 



C 



A 



C 



or 



ri cot — - rs tan — + 2^/^^ = rA cot — - ta tan — 

A C „ , A C 

n cot — - rs tan — - 2^rir3 = r^ cot -^ - ^4 tan — , 



and we obtain from the third disjunction that 



A 



B 



or 



A 



B 



Define u, v, w, x, y, z by 



(25) 
(26) 

(27) 
(28) 



ri cot — 


- r2 tan — + 2y/rir2 = rA cot — 


- rA tan — 


(29) 


A 
n cot — - 


B A 

- r2 tan — - 2^1-1 ^2 = rA cot — - 


B 

-rAtan—. 


(30) 


, m, n by 








A 
l = cot-, 


B 

m = tan — , 


C 
n = tan — . 


(31) 



\/r2r3 

TA 


if (ESI) holds, 


\/r2r3 

TA 


if dMl) holds. 


y/rir3 
rA 


if (127| holds. 


^/rirs 
rA 


if (EHl) holds. 


\/rir2 
rA 


if ^ holds, 


\/rir2 


if (I3()l) holds. 



rA 

ri r2 

X = —, y ^ —, 
rA rA 



rA' 



Then we have 



(32) 



' my + nz ~ 2u ~ fh + n^ 
Ix — nz + 2v = I — n, 
Ix — rhy + 2w = I — rh, 
xy = w'^, 
xz = v'^, 
yz — u . 
In Case 3, we have 

^BD2-DiC + D2D3 = BDa + Da C, 
-AEi + CE3 + E1E3 = AEa - CEa or - AEi + CE3 - E1E3 = AEa - CEa, 
-AFi + BF2 + F1F2 = AFa - BFa or - AFi + BF2 - F1F2 = AFa - BFa- 

By expressing the lengths by the radii and the angle sizes, we obtain from the 
first equation that 



B C „ , B C 

- ^2 tan — - ra tan — + 2^r2r3 = ta tan — + rA tan —, 

we obtain form the second conjunction that 

A C „ , A C 

-n cot — + ra tan — + 2y/nr3 = ta cot -^ - 14 tan — 



or 



A 



C 



A 



C 



-ri cot — + rs tan — - 2^/rIri' = rA cot "^ - 14 tan — 
and we obtain from the third conjunction that 



A 



B 



, A B 

-n cot — + r2 tan — + 2i/rir2 = r^ cot -^ - 14 tan — 



or 



B 



B 



-ri cot — + r2 tan — - 2^/nr2 = rA cot "^ - 14 tan — . 
Define I, fh, n hy (|5T|) . Define u, u, w, x, y, z by 



u = — 


V'^2r3 

rA 




14 


if dMl) holds, 




I 14 


if ([H5l) holds. 




y/rir2 
1 14 


if (|36|) holds. 




Vrir2 

1~A 


if daZl) holds. 




—, y = 

14 


r2 _ 
rA' 


H 

rA 



Then we have the same system of equations as ([5^ . 



(33) 

(34) 
(35) 

(36) 

(37) 



10 



For any triangle ABC, ii I, fh, n are defined by pip, then Ifhn — I — rh — n 
holds. On the other hand, if positive reals I, fh, n satisfy Iffih = I — fh — h, then 
there exists a triangle ABC that satisfies ([31]) . Thus, Case 2 and Case 3 can be 
unified and reduced into solving the system of equations ([5^ for u, v, w, x, y, 
z with positive real parameters /, fh, h under the restriction Ifhh, = I — fh — h. 

As we will show in Appendix |^ the system of equations has the following 8 
solutions. 



^/P 



1-/-1 



\fmF 



1 + m- 1 



\/v?- + 1 + n — 1 



v/^ 


+ 1- 


-Vm^ 


+ 1- 


-Vn^ 


+ 1 + 1- 


fh — h 


+ 1 










21 




-1 + 


fh + h 




v/^ 


+ 1- 


-Vm^ 


+ 1 + Vfi^ 


+ 1- 


-1 










2fh 




-1 + 


fh + h 




v/^ 


+ 1 + Vm2 


+ 1- 


-^n^ 


+ 1- 


-1 



(38) 



2n 



VP 



l-l-l 



\Jvn?- + 1 — 771+1 



\Jv? + 1 — 71+1 



\fW+\ + \Jfh^ + 1 + Vfi^ + 1 + ? -7Ti-n + 1 



(39) 



V = 



21 


- ^ + 777 + 77 - 




Vl^ + l + Vfh^ + 1 - Vr72 + 1 - 


-1 


2?7l 


- / + 777 + 77 - 




%//2 + 1 - V?7l2 + 1 + Vf72 + 1 - 


-1 



277 



Vl2+T-; + i 



V7772+T — 777 — 1 



Vn2 + 1 + 77 + 1 



\/P + 1 + Vm^TT- \Jh^ + 1 + / -777-77-1 



(40) 



y = 



21 




yjr- + \ + ^m? + 1 + Vr72 + 1 - 


- ; + 777 + n + 1 


2777 




V^2 + 1 - ^/fh'^ + 1 - V"2 + 1 - 


-/ + m + 77 + 1 



277 



11 



VpTT-1 + 1 



w 



Vm2 


+ 1 + 


m+ 1 




2 


1 


\ffi 


2 + 1- 


-n- 1 



y/P + 1- Vm^ + T+ y/n'^ + 1 + 1 - rh ~ n - 1 



21 




Vl^ + 1- Vm^ + 1 - Vn^ + 1 - 


-l + rh + n + 1 


2m 




V/2 + 1 + Vm2 + 1 + Vri2 + 1 - 


-l + fh + n + 1 



2n 



Vl^^l + l-1 



(41) 



y/im? + 1 — 771—1 



V'T-2 + 1 — n — 1 



V?2 


+ 1- 


- y/m^ 


+ 1- 


-Vri2 


+ 1-1 + 


m 


+ n 


+ 1 










21 




rh 


— h 




V;2 


+ 1- 


- yjfh^ 


+ 1 + Vn2 


+ 1 + 1- 


-1 










2m 




fh 


— n 




V;2 


+ 1 + V?Tl2 


+ 1- 


-Vn2 


+ 1 + 1- 


-1 



(42) 



2n 



y/P + l + l-1 



\Jm^ + 1 + m+\ 



\/n2 + 1 + n + 1 



VP 


+ 1 + 


y/rn? 


+ 1 + Vfi2 


+ 1-; + 


m 


+ n 


+ 1 








2? 




m 


— n 




V;2 


+ 1 + \/?Tl2 


+ 1 - Vn2 


+ 1 + ?- 


-1 








2m 




m 


— n 




V;2 


+ 1- 


y/m^ 


+ 1 + Vfi2 


+ 1 + ;- 


-1 



(43) 



2n 



12 



VpTT + 1 + 1 



VTO^+T +771—1 



W 



\Jn^ + 1 — n + 1 



V^2 + 1 + ^^2 + 1 - ^n^ + l-; + m + n- 1 



2? 


- m - 




V^2 + 1 + \/m2 + 1 + Vri2 + 1 + ; - 


-n + 1 


2m 


- m - 




V^2 + 1 - V7Tl2 + 1 - Vri2 + 1 + / - 


-fi + l 



2n 



VPTT + / + 1 



V"?+T— 771 + 1 



VtT^+T +77—1 



X = 



\/P + 1- ^Jm? + 1 + V"^ + 1-^ + 771 + 77- 1 



2? 


- 777 - 




V;2 + 1 - V?772 + 1 - V7l2 + 1 + ; - 


-77+1 


2777 


- 777 - 




V^2 + 1 + \/7772 + 1 + V772 + 1 + ; - 


-77+1 



277 



Define aA^f^A^lA G (0,+oo) and <ja by 



sinh (XA 



sinh /3a 



(TA = 



s — a 



sinh 7A 



s-6 



The fourth equation in ((38|) corresponds to a value of ri as follows. 

TA (\/l2~+T - ^777^ + 1 - y/n^TT + i-fn-n + i) 



ri 



21 



Since 



rA 



m 



s{s — b){s — c) 
s — a 
s — c s — a 



s — b s — a 



rA 



rA 



re 



rA 



rB 



(44) 



(45) 



13 



it holds that 



TA {VP + 1 - Vm'^ + 1 - Vn'^ + 1 + 1 -fh-n + l) 



21 



rl bcs ab{s — b) ac(s — c) s(s — b){s — c) \ 

= - (cosh OA cosh /3a cosh 7^ — smh a^ cosh f^A smh 74 

— sinha^ sinh/^^ cosh7yi + cosha^ sinh /3a sinh7yi + 1) 

_ rA(cosh(^A + 7A - aA) + 1) 

2 
= rA sinh^(cr4 - cxa). 

By making similar calculations on every last three equations in (l38l) . (l39)) . 
(|40)) . (|4T|) . (|42|) . (I43l) . (|44)) and (gS]), we obtain the following respective solutions 
in Cases 2 and 3. 

ri = rcosh^ (cTA - a^), 

^2 = re sinh^ (cTA ~ Pa), (46) 

rs = r_Bsinh^(cr4 -7^1). 

ri = r cosh cta , 

r2 = re sinh^(cr4 - 7^), (47) 

rs = rBsinh^((TA -/3a)- 

ri = rcosh^ (cTA -7a), 

J'2 = re sinh cta, (48) 

rs = rs sinh^ (o-A - "a)- 

ri = rcosh^(crA - /?a), 

r2 = re sinh^(crA - ua), (49) 

^3 = rs sinh cta • 

ri = rsinh^ (cTA -"a), 

r2 = re cosh^(crA - Pa), (50) 

ra = rs cosh^ (cta -7a)- 

ri = r sinh cta, 

r2 = re cosh^(CTA - 7a), (51) 

ra = rBCOsh^(CTA - Pa)- 

ri = rsinh^(CTA -7a), 

r2 == reCOsh^CTA, (52) 

rs = r_B cosh^(CTA - oa). 

ri = rsinh^(CTA -Pa), 

7-2 = re cosh^(CTA - ua), (53) 

ra = rs cosh^CTA- 
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4.3 Cases 4 & 5 

Define aB,l3B-,lB G (0, +oo) and as by 



h 

-, oiiiii i^tl — Ti °"^ii W — Ti 

— S — S — 

as + Pb + IB 



• V,2 S-C --20 • - 2 S-a 

smn as — Ti smh ps = -r, smn 73 = 



2 
Analogously to l4.2[ we obtain the following solutions in Cases 4 and 5. 

ri = re sinli^(crB - as), 

r2=rcoiih^((7B- Pb), (54) 

rs = rA sinli^(crB -7b)- 

ri = re sinh^CTs, 

ra = rcosh^(crB - 7b), (55) 

ra = rA sinli^(crB - /3b). 

n = re sinh^(crB - 7b), 

r2 = r cosh^ (Tb, (56) 

rs = rA sinh^ (ctb - ob)- 

ri = rp sinh^ ((Tb - /3b), 

r2 = r cosh (ctb - "b), (57) 

ra == rAsinh^CTB- 

ri = re cosh^ (ctb - ^b), 

r2=rsinh2(crB-/3B), (58) 

ra = rA cosh^(crB - 7b)- 

ri = re cosh ctb, 

r2 = rsinh^ (ctb - 7b), (59) 

r-i = rB cosh^(crB - /3b)- 

ri =: re cosh^(crB - 7b), 

r2 = r sinh ctb, (60) 

ra = rA cosh^(CTB - ctB)- 

»'i = re cosh^(crB - /3b), 

r2 = rsinh^(CTB - Ob), (61) 

»'3 = "TA cosh^ Ctb ■ 

4.4 Cases 6 & 7 

Define ae,/3e,7c G (0,+oo) and ac by 

. , 2 s-h ■ .-i a *~" -1,2 c 

smh OLc — , smh pc — , smh jc = , 

s — c s — c s — c 

ac + /3c +7C 
c^c = 7i ■ 
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Analogously to 14. 2[ we obtain the following solutions in Cases 6 and 7. 

n == vb sinh^((Tc - "c), 
r2 = rAsinli^(crc' - /^c), 
, 7-3 = rcosh^{ac - 7c)- 
ri = vb sinh^ ac, 
r2 = r/isinli^(crc -7c), 
ra = r cosh^ (tTc - /3c)- 
ri = rs sinli^(crc' - 7c), 
r2 = rAsinh^crc, 
r3 = rcosh {ac — Q^c)- 
ri = rB sinli^(crc' - l^c), 
r2 = rA sinh^(crc - ac). 



ra 


= rcosh ac- 




ri 


= rBCOsh^(CTc 


-ac). 


»'2 


= r^ cosh^ ((Tc' 


~/3c). 


?-3 


= rsinh^ (cTc — 


7c)- 


?-i 


= rBCOsh^crc, 




»'2 


= rA cosh^(crc 


-7c), 


^-3 


= rsinh^ (cTc - 


/3c). 


7^1 


= vb cosh ((Tc' 


-7c), 



r-i = Mcosh ac, 
ra = rsinh^(crc' - /3c)- 
ri = vb cosh^(crc - Pc) 
^2 = rA cosh^(CTc - ac) 
ra = r sinh ac- 



(62) 



(63) 



(64) 



(65) 



(66) 



(67) 



(68) 



(69) 



5 Conclusion 

Theorem 1. For any triangle, there exist 32 triplets of circles such that each 
circle is tangent to the other two circles and to two of the sides of the reference 
triangle or their extensions. The radii can he expressed by p7|) - (|24l) . (|46l) - (|53l) . 



Figures [5HM1 illustrate the 32 triplets of circles for a triangle ABC such that 
A = 45°, B = 54°, C = 81°. 



A Solutions of the systems of equations 

In this appendix, we will solve some systems of equations by computing Grobner 
bases. Although it is difficult to compute the Grobner bases by hand, any 
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computer algebra system that can compute Grobner bases should work. 

Proposition 1. The system of equations ^ for the variables u, v, w, x, y, z 
with the positive real parameters I, m, n under the restriction Imn = I + m, + n. 
has 8 solutions (|9))- p6|) . 

Proof. Counting Z, m, n among the variables in addition to u, v, w, x, y, z, we 
compute the reduced Grobner basis of {my + nz + 2u — m — n, lx + nz + 2v — l — n, 
Ix + my + 2w — I — m, xy — w^ , xz — v"^ , yz — u^ , Imn — I — m — n} with the 
degree reverse lexicographical ordering x)^y)~z)~u)~v)^w)^l)^m)^n. 
The reduced Grobner basis consists of 67 polynomials including 

/i = {2u^ + 2lu - 2-u - l){2u^ + 2lu + 2u + I), 

/2 = (2w^ + 2mv -2v~ m){2v'^ + 2m,v + 2u + m), 

/s = (2w^ + 2nw - 2w - n)(2w^ + 2nw + 2w + n), 

fi — Ix — u + V + w ~ I, 

fz — my + u— v + w — m,, 

fe = nz + u + V — w — n. 

By solving {/i = 0, /2 = 0, /a = 0, /4 = 0, /a = 0, /e = 0} for u, v, w, x, 
y, z, we obtain 128 solutions. Note that Imn = I + m + n is equivalent to 

I = ^^^±^. (70) 

mn — 1 

By assigning each of the 128 solutions together with ([701 to {my+nz+2u—m,—n, 
lx + nz + 2v — l — n, lx + my + 2w — l — m, xy — w'^, xz — v^, j/z — m^} and then 
picking out the solutions such that the assignment makes all of the polynomials 
equal 0, we still have 8 solutions (P|)- P^ . which are the solutions of (|S]). D 

Proposition 2. The system of equations p2p for the variables u, v, w, x, y, z 
with the positive real parameters I, m, n under the restriction Imn = I — rh — n 
has 8 solutions (|55|) - (|i5|) . 

Proof. Counting Z, ffi, n among the variables in addition to u, v, w, x, y, z, we 
compute the reduced Grobner basis oi{my + nz — 2u — rh—n, lx — nz + 2v — l + n, 
Ix — ffiy + 2w — I + m, xy — w'^, xz — v"^, yz — u^ , Imn — I -\- fh + n} with the 
degree reverse lexicographical ordering x>y>z>~u>~v>w>~l>fh>~n. 
The reduced Grobner basis consists of 67 polynomials including 

/i = (2u^ + 2lu - 2m - l){2u^ + 2lu + 2u + I), 

/2 = (2w^ - 2mv -2v^ m)(2v^ - 2mv + 2v-m), 

/g := (2^2 _ 2nw ~2uj + n)(2w^ - 2nw + 2w - n), 

fi ~ Ix — U + V + w — I, 

/5 — ffiy — u + v — w ^ rh, 

fe = hz — u — V + w — h. 

By solving {/i = 0, /2 = 0, fs =0, /4 = 0, /s =0, /e = 0} for u, v, w, x, 
y, z, we obtain 128 solutions. Note that Ifrih — I — m — h is equivalent to 

1 = -^^. (71) 

mn — 1 
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By assigning each of the 128 solutions together with dTTI) to {fhy+nz~'2u~'fh—n, 
Ix — hz + 2v — I + h, Ix — ifiy + 2w — I + ffi, xy — w^, xz — v'^, yz — m^} and then 
picking out the solutions such that the assignment makes aU of the polynomials 
equal 0, we stiU have 8 solutions ((55)) ~(P5 |) . which are the solutions of ([5^ . D 
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ri = TA sin^(cr — a), 

Figure 3: { r2 = tb sin^(CT - P), 

rs = rcsm^{a - 7). 




ri = TA sm cr, 

,2/ 



Figure 4: •^ ^2 = r^ sin (cr — 7), 
rs = rcsin^(CT - /3). 
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ri = rASm'^{a- 7), 
Figure 5: { r2 — tb sin^ a, 

ra = rp sin (cr — a). 
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ri — TA sm {G 

,2/ 



Figure 6: \ r-i^ tb sin [a — a), 



ra — re sm cr. 
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fi = TA cos^(cr — a), 

Figure 7: { r2 — tb cos^(cr — j3), 

ry, — re cos^(cr — 7). 
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ri — rA cos (T, 

„2 



Figure 8: \ r2 — tb cos^[a — 7), 
ra = rcCOS^(cr- /3). 
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' ri = rA cos^(ct — 7), 



Figure 9: < 



?'2 — tb cos ct, 
. rs, = rcCOS^{a - a). 
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ri =rACOs^(cr-/3), 
Figure 10: i r2 ~ tb cos^(cr — a), 



rs 



re cos a. 




ri — rcosh {aA — ua), 
Figure 11: ^ r2 = rcsinh^(cr4 - I3a), 
rs = rBsinh^{aA -Ja)- 
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ri = r cosh cta , 
Figure 12: { r2 ^ re sinh^(cr4 - 7^), 
ra = rBsinh^(cr4 - /3a). 



26 




ri = rcosh^ {aA -Ja), 
Figure 13: ^ r2 = re sinh aA, 

rs = rBsinh^(cr4 - a^). 
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ri =rcosh^(cr4 -Pa), 
Figure 14: { r2 — rcsinh^(cr4 — a^i), 
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ri = rsinh^(cr4 - q;a), 
Figure 15: ■( r2 = re cosh^(crA - /3a), 
^•3 = ?'s cosh^((TA - 7a). 
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ri = r sinh (ja , 
Figure 16: { r2 = re cosh^(o>i - 7^), 
ra = rs cosh^(cr4 - 13a). 
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ri = rsinh^(cr4 -7^), 
Figure 17: { r2 = r^ cosh ua, 

rs = rBCOsh^(cr4 - a^). 
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ri = rsinh^(CTA ~ Pa), 
Figure 18: ^ r2 — re cosh^{aA — a_A) 
r3 = rscosh^fJA. 
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ri = re sinh^(crB - as), 
Figure 19: { r2 — rcosh^[aB — /?b), 
ra = rA sinh^(crB -7b)- 
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ri = re sinh ub, 
Figure 20: ^ 7'2 — rcosh^ (ctb — 7s), 
ra = rAsinh^(crB - Pb)- 
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n = re sinh^ {as -1b), 
Figure 21: •^ ^2 = rcosh gb, 

rs = TA sinh^(CT_B - as)- 
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-- rcsinh (ub ~ Pb), 
Figure 22: { r2 ^ rcosh^{aB ~ as), 
rs = rAsinh^ as- 
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ri = re cosh^(o-B -as), 
Figure 23: { r2 = r sinh^ (ctb — /3b), 
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fi = re cosh as, 
Figure 24: { r2 ^ rsinh {as — 7b), 
r-a = rs cosh^ (ctb - Pb)- 
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ri = re cosh^{(7B -7b), 



Figure 25: < r2 — rsinh as, 



ra = rA cosh^(CTB - as). 
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ri = re cosh^(crB - (3b), 
Figure 26: { r2 — rsmh^{aB — as), 
: rA cosh aB ■ 
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Figure 27: { r2 ^ ta sinh^ (ctc - Pc), 
rg = rcosh^(crc - Ic)- 
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>! = rs sinh <Jc, 
Figure 28: { r2 ^ ta sinh^ (ctc ~lc), 
r3 = rcosh^(crc -/3c)- 



42 




ri = rB sinh^ {ac - 7c), 
Figure 29: { r2 — r^ sinh ac, 

r^ = rcosh [ac — etc)- 
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»'i = rs sinh^(crc - /3c), 
Figure 30: ^ ^2 = r^ sinh^(CTc — etc), 
r^ = r cosh ac- 



44 




Figure 31: { 7-2 = rA cosh^(crc' - Pc), 
r-s = rsinh^(CTc - Ic)- 
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fi = tb cosh (7(7, 



Figure 32: 



2 
,2/ 



r2 = rACOsh {ac -Ic), 
ra = rsinh^(CTc' - /?c)- 
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ri = rB cosh^(o-c - 7c), 



Figure 33: 



r2 



rA cosh uc 



ra = rsinh^(f7c' - Pc)- 
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Figure 34: ■^ r2 ~ ta cosh (0-17 -^ cec), 
ra = rsinh ac- 
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